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PLANE  GRILLAGE  SUITE  OF  PROGRAMS 


• 10.^  C W J Simons/  <r„ 

INTRODUCTION  

1 .  A suite  of  plane  grillage  programs  has  been  implemented  on  three  computers 

a.  UNI VAC  1108 


The  features  available  on  each  of  the  computers  are  shown  in  the  table  below 


Computer 


UNIVAC  1108 


Stress  Analysis  (in-core) 
Stress  Analysis  (out  of  core) 
Lowest  Buckling  Mode 
General  Buckling  Modes 
Natural  Frequencies 
Data  Generation 
Data  Plot 


2.  On  all  of  the  computers  full  use  is  made  of  symmetry  to  enable  problem 
size  to  be  reduced. 

3.  The  main  text  of  the  report  is  divided  into  two  sections. 

4.  The  first  section  of  the  report  gives  a brief  summary  of  the  methods  used 
for  the  various  solutions  required,  ie 

a.  Stress  Analysis 

b.  Buckling 

c.  Natural  Frequencies 

The  second  deals  with  the  use  of  the  program.  This  gives  details  of  the  ideal- 
isation of  the  structure  considered,  preparation  of  the  data  for  each  of  the 
facilities  available,  a description  of  the  control  cards  that  are  required  for 
each  of  the  computers  and  details  of  the  output  produced  for  each  of  the 
facilities.  The  detailed  theory  for  the  various  analyses  is  contained  in 
Appendices  A,  B and  C. 

5.  This  report  supersedes  R48l  by  S B Kendrick  and  J L MoKeeman,  which 
describes  an  earlier  version  of  the  program.  It  was  originally  developed  on 
a Pegasus  computer  programmed  in  autcode. 

6.  The  program  is  NCRE  number  N21 , 


I 


GENERAL  METHOD 

7.  A plane  grillage  can  be  idealised  into  two  basic  types  of  finite  beam 

elements.  These  are:- 

k-  » 

i.  Double  - ended  elements  - characterised  by  the  need  to  consider  the 
displacements,  or  forces,  at  both  ends. 

ii.  Single  - ended  elements  - characterised  by  the  need  to  consider  the 
displacements  or  forces  at  only  one  end  (normally  the  second  end  is  at 
a boundary  and  30  has  prescribed  deflections  and  stresses).  Where  two 
or  more  of  these  elements  meet,  a reference  point  is  chosen,  called  a 
hub  position.  Normally  the  finite  size  of  the  beam  cross-section  is 
ignored,  and  in  this  case  the  deflections  and  rotations  of  the  beam 
are  equal  to  those  at  the  hub.  The  program,  though  also  allows  for 
the  general  case  when  the  beams  are  separated  by  a rigid  area,  (see 
figure  1 ) or  when  one  or  more  of  the  beams  axes  do  not  pass  through 
the  hub. 

8.  A stiffness  matrix  for  each  beam  is  formed  using  Timoshenko  beam  theory  (l)* 
allowing  for  shear  deflections.  Torsion  is  also  included.  The  formation  of  this 
stiffness  matrix  is  detailed  in  Appendix  A. 

9.  The  stiffness  matrices  for  the  elements  are  then  transformed  to  the  global 
coordinate  system  and  combined.  See  Appendix  A.  The  resulting  matrix  is  used 
(and  modified  if  necessary)  by  each  feature  of  the  program 

a.  Stress  Analysis  (with  or  without  in-plane  end  loads  coupled  to  initial 

deformations) 

The  set  of  equations  built  up  and  solved  by  the  computer  is 

M ■ tAj  W * W 

[a]  ” overall  structure  stiffness  matrix 

|g  , - loading  matrix  - can  include 

i.  distributed  pressure  loads 

ii.  point  loads 

iii.  initial  deformations  plus  in-plane  end  loads 

^ - force  vector  - externally  applied  moments  and  shear  force 

5 . - displacement  vector  - rotations  and  normal  displacement 

This  equation  is  solved  for  displacements,  and  stresses  are  calculated 
(see  Appendix  A) 

b.  Lowest  Buckling  and  General  Buckling  Modes 


[*]{»}  - \ [»,]  {s} 


[B  ] - aftitrix  representing  the  relative  distribution  of  axial  forces 
' in  the  element,  (see  Appendix  B) 

A - is  the  axial  load  factor  at  which  the  grillage  buckles 


*(  ) = References  on  Rige  16 


w 
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Natural  Frequencies 


reduced  stiffness  matrix  after  G-uyan  reduction  (5)  has  been 
applied  to  [A]  (see  Appendix  C) 

lumped  mass  diagonal  matrix 

natural  frequency  in  radians/sec 


USE  OF  PROGRAM 


Idealising  the  Structure 

10.  The  grillage  must  be  idealised  into  a number  of  structural  elements 
connected  at  hubs. 

a.  Each  hub  is  given  a number  starting  at  1 and  numbered  consecutively, 

b.  Double  ended  elements  are  numbered  at  each  end  using  consecutive 
numbers . 

c.  Each  single  ended  element  is  given  a label  number,  which  must  be 
different  from  anv  other  label  number  used  for  an  element. 

11.  Different  types  of  element  are  available  as  shown  below. 

Type  No 

Double  ended  1 General  beam  between  2 hubs 

2 Beam  with  symmetry  about  one  end 

3 Beam  with  anti  symmetry  about  one  end 

Single  ended  11  Beam  elastically  restrained  at  one  end 

(iL  = -a  V , T = -b  *'  J S,  = -o  * ) 


Beam  simply  supported  at  one  end 


Beam  clamped  at  one  end 


Beam  symmetrical  about  mid-point 

Beam  anti  symmetrical  about  mid-point 

A pillar  support  restraining  transverse 
deflections  but  producing  no  moments  at  the 
point  of  support 


NOTES 


Type  11  - The  values  of  a^ , b^  and  c^  should  be  obtained  from  an 

approximate  estimate  of  the  flexibility  of  the  supporting  structure 

Types  2 and  3 - These  elements  are  included  to  increase  the  size 
of  the  grillage  that  can  be  analysed,  so  that  only  a half  or  a 
quarter  of  the  structure  need  be  considered.  The  loading  must  also 
be  symmetrical  about  any  axis  of  symmetry  and  antisymmetrical  about 
any  axis  of  antisymmetry.  In  using  these  elements,  it  is  only 
permissible  for  them  to  have  one  image  element.  Every  element  of 
the  complete  structure  must  be  accounted  for  in  the  idealisation 
either  as  an  element  or  as  the  image  of  an  element.  For  examples 
of  permissible  idealisations  see  figure  3. 


UNLIMITE 


o.  Type  29  - The  element  produces  a normal  force  to  the  grillage  but  no 
moments.  If  the  parameters  for  this  element  are  set  at  zero  the 
program  assumes  a rigid  simple  support. 

Data  Preparation 

12.  A great  deal  of  the  input  is  common  to  all  the  programs  but  because  of 

minor  variations  individual  specifications  are  given.  The  notation  given  below 

is  common  to  all  the  programs.  Special  notation  for  an  individual  program  is 

given  with  the  notation  for  that  particular  program. 

NT  parameter  to  determine  the  version  of  the  program 

NT  s 1 stress  analysis 

=2:-  lowest  buckling  mode 

= 3;~  general  buckling  modes 
4:-  natural  frequencies 

NBEAM  total  number  of  elements  in  the  grillage 

JOINT  number  of  hubs 

Ivp  Kp>  0 prints  element  stiffness  matrices 

= 0 suppress  such  printing 
N1  label  number  of  element  (at  end  0) 

N2  element  type  number 

HI  hub  number  to  which  end  0 of  element  is  attached 

He  if  double  ended  element,  hub  number  to  which  end  B of  element 

is  attached,  if  single  ended  element  He  = 0 

A effective  shear  area  of  the  element 

I moment  of  inertia  of  cross  section  of  element 

Z distance  of  the  outermost  fibre  of  the  element  cross  section 

from  the  neutral  axis  (positive  if  measured  from  the  flange, 

negative  if  measured  from  the  faying  flange) 

b,  f width  of  faying  flange  and  flange  respectively. 

d depth  of  frame  web 

hb,  hf,  hw  thickness  of  faying  flange,  flange  and  frame  web 
G angle  element  makes  with  the  x axis 

J torsional  constant  (twist  per  unit  length  = Torque/G'J  where 

& is  the  shear  modulus 

0 the  angle  in  degrees  between  element  and  image 

E Young's  modulus 

G Shear  modulus 

B length  of  element  (the  element  is  considered  to  be  labelled  0 

at  one  end  and  B at  the  other  end) 

XI,  Y1  X,  Y coordinates  of  end  0 of  element  relative  to  hub  at  the 

same  end 

X2,  Y2  X,  Y coordinates  of  end  B of  element  relative  to  hub  at 

same  end 

a,  b,  c specified  for  elastically  restrained  elements  (type  11 ) 

Mg  = aVg,  Tg  = -bfg,  Sg  = -cWg 


hb,  hf,  hw 


XI,  Y1 


X2,  Y2 


a,  b,  c 
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A1 

B1 

Vx,  Vy 
W 

Mx,  My 
S 

SIGMAO, 

SIGMAB 

SIGMAM 


NB 


Cross  sectional  area  of  support,  type  29  only 
length  of  support,  type  29  only 
rotations  about  x and  y axis 

displacement  perpendicular  to  plane  of  grillage 
moments  about  x and  y axis 
force  perpendicular  to  plane  of  grillage 

stress  at  end  0 (or  end  B)  of  the  element 

stress  at  a maximum  in  bending  moment  distribution, 
if  any.  (if  there  are  two  maxima  the  larger  of  the 
two  stresses  is  taken) 

SIGMAP  stress  at  point  of  application  of  concentrated  load 

if  any 

The  values  for  the  parameters  can  be  in  any  units  but  they  must  be 
consistent  throughout, 

i.  Stress  Analysis  (with  or  without  initial  deformations) 


a. 


b. 


Format 

NT,  NBEAM,  JOINT,  KP 

1*110 

NLOAD 

Ho 

For  each  element 

(i) 

N1,  N2,  HI,  HC,  NX 

4110,  FI  0.0 

(ii) 

A,  I,  Z , if  Al>0 

3F10.0 

b,  hb,  d,  hw,  f,  lif  if  N1<0 

6fi  0.0 

(iii) 

«,  J,  E,  G 

5F10.0 

(iv) 

B,  XI,  Y1,  X2,  Y2 

5F10.0 

(t) 

p,  w,  m,  for  each  loading 
condition 

3F10.0 

(vi) 

One  set  for  each  loading  condition 
VO,  WO,  VB,  WB  If  NX  / 0 

4F10.0 

(vii) 

For  type  11  elements  only 
A,  B,  C 

3F10.0 

NOTATION  NLOAD 
NX 


m 


VO,  VB 
WO,  WB 


number  of  loading  conditions 

Axial  force  in  the  element,  tension  positive 

Uniform  load  applied  to  element,  positive  downwards 

Concentrated  load  applied  to  element,  positive 
downwards 

A concentrated  load  acts  at  a distance  mB  from 
end  0 of  the  element 

initial  out-of-plane  slope  at  end  0 (and  end  B) 
of  element  (in  degrees) 

initial  out  of  plane  deflection  at  end  0 (and  end  B) 
of  element 
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NOTES  NT  = 1 for  stress  analysis 


If  an  element  has  the  same  loading  data  as  the  one  previous  to  it, 
set  N2  negative  and  specify  up  to  b (iv)  only 

An  example  illustrating  t!:i  s facility  is  given  in  Appendix  D,  numerical 
example  (i),  figure  7. 


ii.  Data  Generation 


This  routine  cuts  down  the  data  preparation  for  the  common  special 
case  of  a rectangular  grillage  of  othogonally  intersecting  beams 
parallel  to  the  coordinate  axes.  The  data  produced  is  suitable  for 
input  to  the  stress  analysis  section  of  the  program. 


There  are 

2 sub-divisions  in  the  data  preparation. 

Section  1 

- Equal  elements  at  uniform  spacing. 

a. 

NX,  NY,  XL,  YL,  NLOAD,  KP 

211/0,  2F1  /■ 
211/0, 

b. 

IB  Cl,  IBC2,  IBC3,  IBC4 

411 0 

c. 

NS1 , NS2  if  IBC2  or  : IBC3  = 3 or  4 only 

211/0 

a. 

For  elements  parallel  to  the  x axis 
A,  I,Z  or  b,  hb,  d,  hw,  f,  hf 

6F1/0./0 

J 

For  elements  parallel  to  the  y axis 

6F1/0./0 

A,  I,  Z or  b,  hb,  d,  hw,  f,  hf 

e. 

SI&X,  SI GY 

2F1J0./0 

f. 

For  each  loading  condition 
ICOND,  PLOAD 

11/0,  F1/0./0 

for  concentrated  loads  only 
(ie  if  ICOND  = -1 ) 

WLOAD,  XX,  YY 

3F1 0.J6 

g. 

Specify  only  if  SI&X  £ 0 and 
SI&Y  £ 0 

AMD(N),  N = 1 , NLOAD 

6F1/0./0 

Section  2 

- Unequal  elements  at  non-uniform  spacing. 

a. 

As  for  Section  1 

b. 

As  for  Section  1 

c. 

As  for  Section  1 

d. 

For  elements  parallel  to  the  x-axis:- 

As  for  Section  1,  but  if  elements  parallel  to  the  x-axis 
have  different  cross-sectional  or  frame  parameters  from 
each  other,  set  A or  B negative  and  specify  A,  I,Z  or 
b,  hb,  d,  hw,  f,  hf  for  each  of  the  remaining  elements 
parallel  to  the  x axis. 

For  elements  parallel  to  the  y-axis :- 

As  for  Section  1,  but  if  elements  parallel  to  the  y-axis 
have  different  cross-sectional  or  frame  parameters  from 
each  other,  set  A or  b negative  and  specify  A,  I,Z  or 
b,  hb,  d,  hw,  f,  hf  for  each  of  the  remaining  elements 
parallel  to  the  y-axis. 


1 
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e 


6F1  f6.f6 


As  (e)  for  Section  1 
As  (f)  for  Section  1 
As  (g)  for  Section  1 


NOTATION 


NX  - number  of  beams  parallel  to  x-axis 
NY  - number  of  beams  parallel  to  y-axis 
+XL  - length  of  beam  parallel  totbe  x-axis 
^YL  - length  of  beam  parallel  to  the  y-axis 
NLOAD  - number  of  loading  conditions 

KP  - parameter  determining  printing  of  stiffness  coefficients 
from  stress  analysis  program  and  listing  of  nodal  and 
element  coordinates  from  the  plotting  program 

= 0 do  not  print  coefficients 

= 1 Print  coefficients 

= 2 print  coefficients  from  stress  analysis  program  and 
list  coordinates  from  plotting  program 

IBC1,  IBC2,  IBC3,  I3C4  - parameters  indicating  the  boundary 
conditions  for  each  of  the  four  sides  of  the  grillage 
symmetry  or  asymmetry  may  only  occur  about  the  upper  and 
right  hand  edge  of  the  grillage 

IBC1  refers  to  the  left  hand  edge  of  the  grillage 

IBC2  refers  to  the  upper  edge  of  the  grillage 

IBC3  refers  to  the  right  hand  edge  of  the  grillage 

IBC4  refers  to  the  bottom  edge  of  the  grillage  (see  figure  4) 

Parameters:-  = 1 for  simple  support 

= 2 for  clamped 

= 3 for  symmetry/asymmetry  about  centre  of  element 
= 4 for  symmetry/asymmetry  about  end  of  element 
NS1  = 1 if  there  is  symmetry  about  the  upper  edge 
= -1  if  there  is  asymmetry  about  the  upper  edge 
= 0 otherwise 

NS2  = 1 if  there  is  symmetry  about  the  right  hand  edge 
= -1  if  there  is  asymmetry  about  the  right  hand  edge 
= 0 otherwise 

ICOND  = 1 for  uniform  load 

= -1  for  concentrated  load 
PLOAD  = uniform  load 

WLOAD  = concentrated  load  applied  to  the  element,  positive 
downwards 


UNLIMITED 


UrfirTMiTED 


XX  = x coordinate  of  concentrated  load 
YY  = y coordinate  of  concentrated  load 

SIGX  = axial  stress  in  beams  parallel  to  the  x axis,  tension  positive 

SIGY  = axial  stress  in  beams  parallel  to  the  y-axis,  tension  positive 

AMD  = out  of  plane  deflection  at  mid-point  of  grillage  (see  notes 
below) 

ELX  (i),  1=1,  N where  N is  the  total  number  of  elements  parallel  to 
the  x-axis 

= lengths  of  all  elements  parallel  to  the  x-axis 

ELY  (i),  1=1,  N where  N is  the  total  number  of  elements  parallel 
to  the  y-axis 

= lengths  of  all  elements  parallel  to  the  y axis 


(i)  the  x and  y coordinates  are  measured  from  the  origin  where 
the  first  hub  is  considered  to  be  situated 


(ii) 


(iii) 


For  Section  2 (d),  (e),  (f)  the  data  for  all  double  ended 
elements  parallel  to  the  x (or  y)  axis  should  be  specified 
first  starting  from  the  beam  parallel  to  the  x asis  at 
y = 0 (or  parallel  to  the  y axis  at  x = 0)  and  specifying 
the  elements  in  order  of  x (or  y)  increasing.  Repeat  for  all 
subsequent  beams  parallel  to  the  x axis  in  the  order  of  y 
increasing  (or  beams  parallel  to  the  y axis  in  the  order  of 
x increasing).  The  data  for  single  ended  elements  parallel 
to  the  x axis  (or  single  ended  elements  parallel  to  the  y 
axis)  should  then  be  specified  for  each  beam  in  the  same 
order  as  for  double  ended  elements. 

The  parameter  AMD  gives  the  out  of  plane  deflection  of  the 
centre  of  the  grillage.  The  first  hub  of  the  grillage  is 
the  origin  of  the  coordinate  system.  The  deflection  at  any 
point  x,  y is  determined  by:- 

WD  = AMD  Sin  Sin  *(y+Y1  ^ 


TT 


XI  is  the  distance  from  the  left  hand  side  of  the  grillage 
to  the  first  hub 

Y1  is  the  distance  from  the  lower  edge  of  the  grillage  to 
the  first  hub 

The  out  of  plane  slope  at  any  point  x,  y is  determined  by 
differentiating  WD  with  respect  to  the  appropriate  direction. 


iii  Lowest  Buckling  Mode 

a.  NT,  NBEAM,  JOINT,  KP 

b.  MAXIT,  ACC 

c.  For  each  element 

(i)  N1 , N2,  HI,  HC,  Nx 

(ii)  If  N1>0,  A,  I,  Z 

If  N1<0,  b,  hb,  d,  hw,  f,  hf 

(iii)  e,  J,  j6,  E,  & 
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Format 

AI10 

110,  FI  0.0 

AI10,  FI  0.0 
3F10.0 
6F10.0 
5F10.0 


5F10.0 


NOTATION 

MAXIT 

ACC 


(iv)  B,  XI,  Y1 , X2,  Y2 

(v)  For  type  11  elements  only 

A,  B,  C 

= maximum  number  of  iterations 
= accuracy  test  for  the  buckling  load 


3F10.0 


Nx  = axial  force  in  the  element,  tension  positive 
NOTES 

(i)  NT  = 2 for  the  lowest  buckling  mode, 

(ii)  Under  in-plane  end  loading  ignoring  buckling,  each  element 
in  the  grillage  will  develop  a longitudinal  force  Nx.  For 
a given  set  of  end  loads  these  forces  can  be  calculated  by 
other  programs,  (eg  NCRE  Program  N22).  If  all  the  end  loads 
are  increased  by  a constant  factor  the  values  of  Nx  will 
all  increase  by  the  same  factor.  For  a given  relative 
distribution  of  values  of  Nx  this  program  will  calculate  the 
actual  values  of  Nx  (having  the  same  relative  distribution 

as  the  specified  Nx)  at  which  the  grillage  will  first  buckle. 
The  associated  buckling  mode  shape  is  also  calculated. 

(iii)  An  example  illustrating  this  facility  is  given  in  Appendix  D, 
numerical  example  (ii),  figure  8. 

iv.  General  Buckling 

Format 


a. 

NT, 

NBEAM, 

JOINT,  KP 

4110 

b. 

GP, 

GQ,  M 

2F10. 0,110 

c. 

For 

each  element 

(i) 

N1 , N2,  HI,  HC,  rix 

4I10,F10.0 

(ii) 

If  N1>0;  A,  I,  Z 

3F10.0 

If  N1< 0;  b,  hb,  d,  hw,  f,  hf 

6F1 0.0 

(iii) 

6,  J,  0,  E,  G 

5F10.0 

(iv) 

B,  XI,  Y1 , X2,  Y2 

5F10.0 

(v) 

For  type  11  elements  only 

A,  B,  C 

3F10.0 

NOTATION 

GP,  GQ  - these  are  the  limits  within  which  buckling  Factors  (xi) 
are  required  (order  immaterial) 


M = number  of  buckling  loads  to  be  calculated  within  the 

limits  GP, GQ 

mode  shapes  are  also  calculated  if  M>0  but  not  if  M<0 

NOTES 

i.  NT  = 3 for  general  buckling  modes 

ii.  An  example  illustrating  this  facility  is  given  in  Appendix  D, 
numerical  example  (iii),  figure  8. 
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Natural  Frequencies 


Format 


a. 

NT,  NBEAM,  JOINT,  KP 

4110 

b. 

GP,  GQ,  M,  EX,  GRAV 

2F10.0 

2F10.0 

\ c* 

For  each  element 

(i)  N1 , N2,  HI,  HC 

4110 

'} 

i 

(ii)  If  N1 >0;  A,  I,  Z 

3F10.0 

1 

If  N1<0;  b,  hb,  d,  hw,  f,  hf 

6F10.0 

i 

(iii)  ©,  J,  0,  E,  G 

5F10.0 

i 

(iv)  B,  XI,  Y1 , X2,  Y2 

5F10.0 

. 

I 

.1 

(v)  jx,  w,  m 

(vi)  For  type  11  elements  only 

3F10.0 

A,  B,  C 

3F10.0 

NOTATION 


GP,  GQ  - Limits  within  which  the  squares  of  the  natural  frequencies  are 
required  (order  immaterial) 

M i.  If  M<0;  up  to  M natural  frequencies  are  calculated,  but  the 

mode  shapes  are  not  calculated, 

ii.  If  M>0;  up  to  M natural  frequencies  and  their  associated 
shapes  are  calculated 

EX  - A scaling  factor,  (if  left  at  zero  a value  of  3 x 10^  is  set 
within  the  program)  Its  value  must  be  of  the  order  of  Young's 
Modulus 

9 

GRAV  - acceleration  of  gravity  eg  386,4  in/sec 
H - Mass/unit  length  of  an  element 

w - Lumped  mass 


Lumped  mass  is  positioned  at  a distance  mB  from  end  0 of  the 
element 


NOTES 


i,  NT  = 4 for  natural  frequencies 

ii.  An  example  illustrating  this  facility  is  given  in  Appendix  D, 
numerical  example  (iv),  figure  9# 

vi.  Stress  Analysis.  Normal  Loading.  Disk  Solution  for  larger  Problems 

a.  NBEAM,  NLOAD,  JOINT,  KP  4110 

b.  (i)  IB  AND,  MAXIB  2110 

(ii)  KBAND(I),  1=1,  IBAND  6110 

c.  For  each  element 

(i)  N1,  N2,  HI,  HC  4110 

(ii)  If  N1>0;  A,  I,  Z 3F10.0 

If  N1<0,  b,  hb,  d,  hw,  f,  hf  6F10.0 

(iii)  «,  J,  ft,  E,  & 5F10.0 
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(iv)  B,  XI,  Y1,  X2,  Y2 

(v)  For  each  loading  condition 
p,  w,  m 

(vi)  For  type  11  elements  only 
A,  B,  C 


5F10.0 


3P10.0 


3F10.0 


NOTATION 


NLOAD  - number  of  loading  conditions 

IBAND  - number  of  bands  into  which  the  grillage  is  divided  (see  below) 
MAXIB  - largest  number  of  hubs  in  a band 
KBAND  - highest  hub  number  in  each  band 

p - uniform  load  applied  to  element,  positive  downwards 
w - concentrated  load  applied  to  element,  positive  downwards 


- a concentrated  load  acts  at  a distance  mB  from  end  0 of  the 
element 


NOTES 


i.  The  equation  solution  routine  uses  banded  matrix  techniques. 

Storage  of  the  equations  on  disk  backing  store  enables  solution 
of  up  to  about  24.0  degrees  of  freedom  (3  degrees  of  freedom 
per  hub) 

ii.  The  hubs  must  be  numbered  in  bands  such  that  a beam  in  the  nth 
band  can  only  link  to  hubs  in  that  band  or  hubs  in  the  (n+1 )th  band. 

iii.  The  program  is  most  efficient  for  smaller  bandwidths.  For  example, 
a rectagular  grillage  would  be  most  efficiently  banded  if  the  hubs 
were  numbered  in  strips  taken  along  the  shorter  of  the  two  sides. 

iv.  If  an  element  has  the  same  loading  data  as  the  one  previous  to  it, 
set  N2  negative  and  specify  as  far  as  c (iv)  only. 

vii.  Plotting  of  a Plane  Grillage 

The  input  for  this  program  is  in  the  form  of  data  used  for 
i.  and  vi.  At  the  end  of  the  data  the  following  parameters 
must  be  inserted. 


a.  + PL 

b.  If  PL<0 

A,  IT 


FI  0.0 


A8,  110 


NOTATION 


PL  - The  length  in  cm.  to  which  the  user  wants  the  grillage  plotted. 
This  length  does  not  include  space  for  notation  on  the  plot  which 
takes  up  to  4 cm.  of  the  usuable  paper  width. 

A - Name  of  the  file  which  will  hold  the  plotting  instructions 

IT  - The  plotting  time  in  minutes 

NOTES 

i.  This  program  is  used  chiefly  for  checking  that  the  data  is  correct. 

13.  There  follows  some  general  comments  on  data  preparation. 

a.  In  order  to  specify  the  minimum  amount  of  data,  use  is  made  of  the 
sign  and  value  of  N2,  the  type  number. 
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i.  N2  = 0 is  specified  for  an  element  which  is  identical  in  all 
respects  to  the  previous  element  (either  double  ended  or 
single  ended).  In  this  case  only  N1,  0,  HI,  HC  are  specified 
for  the  element. 

ii.  N2  positive  is  specified  for  an  element  whose  data  is  in  the 
normal  form. 

iii.  N2  negative  is  specified  for  an  element  which  has  some  of  its 
data  the  same  as  the  previous  element.  The  loading  conditions 
must  all  be  the  same  as  the  previous  element  if  this  facility 
is  to  be  used.  It  is  necessary  to  specif^  only  up  to  and 
including  the  last  value  which  is  different  from  the  previous 
element  with  the  exception  that  A,  I,  Z or  b,  hb,  d,  hw,  f,  hf 
are  considered  as  groups  and  if  any  of  the  group  is  different 
from  that  of  the  previous  element  then  the  whole  group  must 

be  specified.  A special  case  of  this  is  to  specif^  only  N1 , 
N2,  HI,  HC  where  N2  is  negative.  This  gives  an  element  which 
is  identical  in  all  respects  to  the  previous  element  but  has 
a different  type  number. 

b.  If  required  certain  standard  values  are  set  in  the  program.  The 
standard  value  of  a parameter  is  used  in  the  program  if  a zero 
value  of  the  parameter  is  specified  in  the  data.  The  values  used 
are:^ 

E = 3 x 107  psi 
& = 1 .154  x 107  psi 
0 = 180° 

JiO 

XI  = X2  = Y1  = Y2  = 0 
A1  = 10  sq  in 
B1  = 1 in 

Control  Cards 

14.  A description  of  the  control  cards  for:- 

a.  IBM  1130 

b.  UNIVAC  1108 

c.  KDF-9  follows 


a.  IBM  1130 

i.  Control  for  stress  analysis,  lowest  buckling  mode,  general  buckling 
modes  or  natural  frequencies 

//bXEQb  KCHItfcWM 

•L0CALKCH1 , CSE1 , CSE2 , CSE1 1 , CSE1 2 , CSE1 4, CSE1 5, CSE29, CSJRM, CBIGM,CNCS0 
. | Data 

ii.  Control  for  data  generation 
//bXEQBDN21D 

:{ 

. < Data 
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b.  UNIVAC  1108 


Control  for  stress  analysis,  lowest  buckling  mode,  general 
buckling  modes  or  natural  frequencies 

<§RUN  Jobname,  Acc  No,  Project  No 

<§&SG,T  T.,  U9,  XXXXXR 

@ASG,T  Z. 

<§M0VE  T.,f> 

(§C0PIN  T.,Z. 


F1.,Fi^S  ///1000 
F2. , FU0  ///'  1000 


@ASG,T 


i.  Control  for  stress  analysis  (in-core  solution),  and  lowest 
buckling  mode 

•JOB 

*XEQ 

•I0DUNIT1  /DI  SC/ 300 
*I0DUNIT2/DISC/3^ 

•I0D0NIT5/READ 

*iodunit6/fkemt 

*STQRAGi//26A00 
•DISCFR0GRAMNCREN21 /NCREUSE1 
•DATA 


•ENDJOB 

ii.  Control  for  stress  analysis  (out-of-core  solution) 
•JOB 
•XEQ 

•I 0DUNIT1 /DISC/300 
•I0DUNIT2/DISC/300 
•I  OD  UNIT  5/READ 
•I0DUNIT6/FRINT 
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•storage/ /26400 

•DISCH10GRAMNCREN21 B//NCREOSE1 
•DATA 


•ENDJOB 

iii.  Control  for  Data  Generation 

•JOB 

*XEQ 

•I ODUNIT 5/READ 

•I0DUNIT6/HHMT 

•DI SC  FR0GRAMNCREN21 D//NCREUSE1 

•DATA 


Data 


•ENDJOB 

iv.  Control  for  Data  Plot 

•JOB 

•JCEQ 

•IODONIT5/READ 
•IODONIT6/PRINT 
•SOBFILEARLLIBRARY 
•DISCPR0GRAMNCRE21 H//NCREUSE1 
•DATA 


•ENDJOB 

Description  of  Output 

15.  1.  Stress  Analysis  (in  Core  and  Out  of  Core 


Displacement  and  forces  - for  each  loading  condition  the  hub  dis- 
placements (Vx,  Vy,  W)  are  printed  preceded  by  their  appropriate 
hub  number.  Next  follow  the  internal  beam  forces  at  hubs  (Mx, 

My,  S)  preceded  by  the  appropriate  hub  and  element  label  numbers. 
For  each  loading  condition  a similar  table  of  displacements  and 
forces  is  given.  For  hubs  on  one  or  more  axes  of  symmetry  or 
antisymmetry  the  hub  forces  and  moments  due  to  the  image  are  also 
included  so  that,  for  example,  moments  about  an  axis  of  symmetry 
appear  xero  (or  very  small)  even  although  the  moments  in  the 
beams  on  each  side  of  the  axis  of  symmetry  are  non-tero.  The 


r F 
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shear  forces  at  such  positions  are  double' tfea ' Vii$es HStr^SS‘ 

single  beams.  If  moments  or  shear  forces  are  required  at 
such  positions  they  are  easily  calculated  from  the  forces  and 
moments  at  the  other  end  of  the  beam  and  the  loads  acting  on 
the  beam.  At  axes  of  antisymmetry  the  situation  is  reversed 
and  moments  are  doubled  whereas  shear  forces  appear  as  zero 
or  very  small. 


b.  Stress  - After  the  tables  of  hub  displacements  and  forces 
for  every  loading  condition  a table  oi  stresses  (tension 
positive)  is  given  in  five  columns  as  follows 

KiiEMSNT  - the  label  at  end  0 of  the  beam 

SIGMA  0 - the  stress  at  the  end  of  the  beam  corresponding 
to  the  label  number 

SIGMA  B - the  stress  at  the  other  end  of  the  beam 

SIGMA  Max  - the  stress  at  a maximum  (dM/dx  = 0)  in  bending 
moment  distribution^if  any 

SIGMA  P - the  stress  at  the  point  of  application  of  the 
concentrated  loadfif  any 

ii.  Lowest  Buckling  Mode 

The  lowest  buckling  load  is  printed,  followed  by  the  associated 
buckling  mode  shape. 

iii.  General  Buckling  Modes 

Up  to  M buckling  loads  are  printed,  each  followed  by  its  associated 
shape. 

iv.  Natural  Frequencies 

Up  to  M natural  frequencies  (rad/sec)  are  printed,  each  followed  by 
their  associated  shapes. 

v.  Data  Generation 

Line  printer  output  of  the  prepared  data  is  produced  for  use  with  the 
stress  analysis  (in  core)  program.  The  IBM  1130  version  also  produces 
a punched  tape  of  the  data. 

vi.  Data  Check 

a.  Line  printer  output:-  If  the  parameter  RP  in  the  initial  data 
is  set  equal  to  2,  the  nodal  coordinates  (X,  Y)  and  the 
element  coordinates  (Xj6,  Yj6,  XB,  YB)  will  be  listed. 

X - X coordinate  of  a hub 

Y - Y coordinate  of  a hub 

X0-  X coordinate  of  end  0 of  an  element 

Y0-  Y coordinate  of  end  0 of  an  element 

XB-  X coordinate  of  end  B of  an  element 

YB-  Y coordinate  of  end  B of  an  element 

b.  Plotter  Output:-  Lines,  denoting  the  beams,  are  drawn  with  hub 
numbers  and  element  type  numbers  at  the  boundaries.  Uniform 
loading  is  shown  by  2 small  'v's  positioned  at  each  end  of  the 
element,  and  a concentrated  load  is  shown  by  one  large  *v'  at 

a distance  mB  from  end  0 of  the  element. 
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CONCLUDING  REMARKS 

16.  Sample  problems  have  been  run  and  timed  on  the  various  computers  and  these 
running  times  can  be  found  in  Appendix  E. 

17.  There  are  different  limitations  for  the  different  computers  on  the  size  of 
problem  that  can  be  analysed.  Restrictions  on  the  size  of  problem  are  set  out 
in  Appendix  F.  Certain  elements  are  not  available  for  certain  facilities  and 
details  of  this  restriction  can  also  be  found  in  Appendix  F. 

18.  Appendix  & gives  a general  flow  chart  for  the  suite  of  programs,  giving 
details  of  stress  analysis,  lowest  buckling  and  general  buckling  loads  and 
modes  and  natural  frequencies. 
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Element  Stiffness  Matrices 

1.  The  equations  to  be  solved  are  as  follows. 


w- 

[A,] 

M 

♦ t»2]  {&c}  * {», 

M 

= [a3] 

jsh 

[ ♦ [»4]  {sc}  . {t 

>2} <2> 

A1 

■ <■«> 

i = 

1,3; 

3 = 1,3 

A2 

s <*«> 

i = 

1,3; 

3 = 4,6 

A3 

- <y 

i = 

1,3; 

3 = 4,6 

\ 

= (b.  .) 

i = 

1,3; 

j = 1,3 

B1 

- <*,t> 

i = 

1,3 

B2 

-<V 

i = 

1,3 

*"i 


are  the  stiffness  matrices  whose  coefficients  are  formed  from 


Timoshenko  beam  theory  (see  l).  The  coefficients  are  as  follows. 
&11  = b^  = El  (4+1 2p)/B(l  +1 2p ) 

*12  * *12  ' 0 

a13  = -b13  = 6EI/B^(l+I2p) 

*21  * b21  “ ° 

*22  * *22  * M/B 

*23  = b23  = ° 

*31  * -t3t  = *13 
a = b =0 
32  32 

a33  = b33  = 12EI/B3  (1+1 2P) 
a14  = b14  = 2EI  (1-6P)/B(1+12P) 

a15  =b15  = ° 
al6  = _bl6  = "al3 
a24  = b24  = ° 
a25  = b25  = ”a22 
a26  - b26  = 0 
a34  = b34  = a13 
a35  = b35  = ° 
a36  = b36  = "*33 

a17  = -(PB2/12)  -WB(l-m)  m[(l-m)  +ffc]/(l+12p) 

= (pB2/l2)  + WB  (l-m)  m [ m+6(3]/(l+12P) 


a^  = -(pB/2)  - W(l-m)  [(l+2m)  (l-m)  +1 2p]/(l  +12 p) 
b?7  = -(pB/2)  - Wm[m(3-2m)  + 12p]/(l+12p) 

where  = El/B^  GA 

NB  The  B matrices  are  the  loading  conditions  on  the  grillage* 
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and  j M l are  the  force  vectors  at  ends  h and  c of  a beam  element.  They 


consist  of  the  bending  force,  the  torsional  force  and  the  shearing  force, 
jg^  jand  jg^  jare  the  displacement  vectors  at  the  ends  h and  c of  a beam 

element. 

Forces  and  displacements  specified  so  far  refer  to  the  local  coordinate 
system  of  each  element.  The  local  forces  and  displacements  must  be  transformed 
to  the  global  coordinate  system.  See  figure  5. 

Equations  (l ) and  (2)  can  be  combined  into 


Befine 


(“4 


2 

\ 


B„ 


as  being  the  force  vector  in  the  global  coordinate  system  and 


w 


as  displacements  in  that  system. 

Define  a transformation  matrix 


T = 


Cos  0 
-Sin  0 
0 


Sin  0 
Cos  0 
0 


0 is  the  angle  that  the  beam  makes  with  the  x-  axis  of  the  global  coordinate 
system. 

i r n r 

6h 
be 


fc] 

[“] 


[ 


TA 


TA„ 


TA, 


TA, 


1 

( 


0 

T 


TB 


1 


TB, 


] 


n-1 


- m 

m 

. *■ 

6xh 

T 

0 

B, 

B2 

6xc 

+ 

0 

T 

- - 

- 

« 

(3) 


The  equations  are  put  in  a form  which  can  be  solved  by  equating  to  zero  the 
sums  of  the  moments  and  forces  at  every  hub.  This  provides  a sufficient  set  of 
equations  to  find  the  unknown  hub  rotations  and  displacements.  ", 

If  the  hub  positions  are  not  at  the  end  of  the  beam  it  is  necessary  to 
transfer  the  origin  to  points  distance  x^y^  and  *^2  from  ends  h an<*  c respectively. 

Define  force  vectors  at  the  offset  hub  as 


and  |mxc| 

Define  displacement  vectors  at  the  offset  hub  as 


and  j6x(^ 


|8x»  and 

Define  a transformation  matrix  for  the  forces  at  the  ends  of  the  beam 


T1  - 


1 

0 

0 


0 

1 

0 


-x 

y 
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APPENDIX  B 
THEORY  FOR  BUCKLINO 
The  Differential  Stiffness  Matrix 
The  equation  to  be  solved  for  buckling  is 

a11  a12  a1 3 al4  a1 5 ai6  Txh 

a21  a22  a23  a24  a25  &26  vyh 

a31  a32  a33  a34  a35  *3*  Wh 

ai  . a!  o a.  7 ai  I ai  c ai  C V 

41  42  43  44  45  4°  xc 

a51  a52  a53  a54  a55  &56  vyc 

_a6l  a62  a63  a64  a65  a66  ^ Wc 


(^ij) 


is  the  differential  stiffness  matrix. 


For  details  of  how  it  is 


formed  see  (6).  The  coefficients  are  as  follows 


: 2NxB/l5 
= Nx/10 
:*NxB/30 
:-Nx/lO 
: Nx/10 
: 6Nx/5B 
: Nx/10 
:-6Nx/5B 

r-NxB/30 
« Nx/10 
: 2NxB/l  5 
:*Nx/l  0 
:-Nx/l  0 


where  Nx  is  the  axial  force  in  the  element,  tension  positive  and 
B is  the  length  of  the  element. 


Errors  in  the  Critical  Load  (see  7) 

a.  Beam  simply  supported  at  either  end:- 

i 

If  only  one  element  is  used  to  model  the  beam  an  error  of  approximately 
20$  is  present  in  the  critical  load.  In  the  case  of  this  program  this 

r j 

problem  is  academic  because  of  the  facility  of  single  ended  elements. 
Two  elements  are  needed  to  model  a simply  supported  beam  and  for  this 
case  ''the  error  is  less  than  1$. 

b.  Beam  clamped  at  one  end  and  free  at  the  other. 

For  one  element  the  error  in  the  critical  load  was  found  to  be  0.73f9? 
(according  to  (7)  the  error  is  0,752$).  For  more  than  one  element 
the  error  was  found  to  be  less  than  1$. 


ii.  Evaluation  of  Lowest  Buckling  Mode 

The  equation  to  be  solved,  is 


[A]  is  the  same  as  the  [A]  matrix  derived  in  Appendix  A. 
[BD]  is  the  differential  stiffness  matrix 
To  solve  the  equation  pre  multiply  throughout  (l ) by  [a] 

([ I]  - [U]  X,)  [6]  = 0 where  [u]  = [a]”1  [bD] 


and[j]  is  the  unit  matrix 
Let  |J  = lA,  then  we  get 

(P  [I]  - [U])  [5]  = 0 (2) 

This  has  now  been  reduced  to  the  standard  eigenvalue  problem. 

The  root  of  interest  of  the  equation  p [I]  - [U]  is  the  dominant  root  since 
this  gives  the  smallest  value  of  Xj • The  dominant  root  can  be  found  by 
repeated  by  pre-multiplication  of  an  arbitary  column  by  [U],  (For  theory  see 

Evaluation  of  General  Buckling  Modes 

Again  we  are  solving 


This  can  be  transformed  into  the  standard  symmetric  eigenvalue  problem 


by  carrying  out  a Choleski  decomposition  of  B into  B = LL  (see  3)  and  setting 
-1  -T  t 

E = L AL  and  x = L 6»  This  standard  problem  is  solved  in  the  program  by  an 
HGW  type  of  method  (see  4). 


Appendix  B3 

UNUMITE 


APPENDIX  C 


1 . 


THEORY  FOR  NATURAL  FREQUENCIES 

The  vibration  problem  can  be  written  in  matrix  form  as  follows. 


W2  [M][6]  = [A]  [6] (1) 

The  [Ajmatrix  is  the  same  as  that  derived  in  Appendix  A 

The  [m]  matrix  is  the  same  as  the  [b]  matrix  derived  in  Appendix  A 
except  that:- 

p = mass/unit  length  of  an  element 
w - lumped  mass 

m = the  lumped  mass  acts  at  a distance  mB  from  end  0 of  the  element 
The  equation  (l)  can  be  written  in  the  following  form. 


where  A^ 


where  8 is  a vector  of  rotations  and  J is  a vector  of  displacements. 
Rotary  inertias  have  been  neglected  because  they  are  very  small  compared 
to  translational  inertias. 


A 9 + A2  W = 0 
^6  = - A^1  A2  W 

also  m 5 + [A4  - A3  A^1  A2]  £ = 0 

Rotations  have  been  eliminated  and  the  final  equation  to  be  solved  is 


(E  - u>2  *)  W = 0 
[E]  W = u)2  [■]  JJ 

This  is  the  Cuyan  reduction  (see  5) 

The  $ vector  now  defines  the  relative  displacements  in  the  vibration  mode 
corresponding  to  the  natural  frequency.  u> 

As  in  the  buckling  case  we  can  factorise  [E]  by  Choleski  Factorisation 
[E]  = [L]  [L]T 

* ([l]-'  [.]  [trT>  <(dT  (»])  = i/»2  <[L]T  [W]  ) 

2 

Therefore  the  matrix  has  eigenvalues  Xi  = t/wi  and  eigenvalues  qi  such 
that  qi  = [l]T  [Wi] 

The  eigenvalue  problem  is  now  solved  using  a standard  HGW  method. 
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APPENDIX  D 
NUMERICAL  EXAMPLES 

1.  Four  numerical  examples  are  shown 

i.  Stress  analysis  for  grillage 

ii.  Lowest  buckling  mode  for  a grillage 

iii.  General  buckling  modes  for  the  same  grillage  as  in  ii. 

iv.  Natural  frequencies  of  a free-free  beam 

2.  Figures  7-9  show  the  idealisation  used  for  each  of  these.  Listing  of  the 
data,  and  the  results  obtained  from  these  examples  are  as  follows:- 


5?  ■ ' < 


LOADING  CONDITION  V 

HOB  VX  v Y M 

1 0.10868756E-06  -0.91556011E-05  0.53  71  1 726E-03 

2 0.60793356E-05  -0. 1 82  S 140  7E-06  0*107230988-02 

3 0*000000006  00  -0.2 1266997E-06  0.12675671E-02 


HUB  ELEMENT  MX  MY  3 


1 

X 

X 

0.13285990E  03 

O.OOOOOOOOE  00 

-0.36665081E  02 

1 

X 

3 

•0.X3286231E  03 

O.OOOOOOOOE  00 

0.39825827E  02 

1 

X 

8 

O.OOOOOOOOE  00 

0.6156581 7£  03 

0.28839228E  02 

X 

9 

o*ooooooooe  00 

-0.61565258E  03 

-0.2500C003E  02 

i 

2 

3 

-0.68636633E  03 

O.OOOOOOOOE  00 

-0.21 138922E  02 

i 

2 

3 

0*888359376  03 

O.OOOOOOOOE  00 

—0.2  7 588 71 3£  02 

j 

2 

7 

0*00000000E  00 

0. 18369988E  08 

0.737 19650E  02 

2 

XO 

0.00000000E  00 

-0.18369929E  08 

— 0.250000C3E  02 

* 

3 

8 

-0*33988165E  03 

O.OOOOOOOOE  00 

-0.328 15298E  02 

3 

8 

0.00000000E  00 

0.17055180E  08 

0.89631073E  02 

3 

XX 

0*000O0000E  00 

-0.17055161E  08 

-0.25000003E  02 

LOADING 

CONDITION  2 

HUB 

VX 

VY 

W 

X 

0.13653952E-07 

— 0.6291 98 38E-08 

0.36913030E— 06 

2 

0.25880018E-07 

-0.27558222E-07 

0.16167896E— 05 

3 

O.OOOOOOOOE  OC 

-0.76121381E-07 

0.29071870E-05 

HUB  ELEMENT 

MX 

MY 

S 

X 

X 

-0. 288765 X8E  01 

O.OOOOOOOOE  00 

—0.2065808  3 E—0 1 

X 

3 

0*288768 75E  01 

O.OOOOOOOOE  00 

-0.133 185 18E-01 

3 

X 

8 

O.OOOOOOOOE  00 

0.56628072E  00 

0.33976608E— OX 

X 

9 

O.OOOOOOOOE  00 

-0.56627857E  00 

O.OOOOOOOOE  OC 

i 

2 

2 

0*1208 X675E  01 

O.OOOOOOOOE  00 

0.20656079E-01 

2 

3 

-O.X208X775E  OX 

O.OOOOOOOOE  00 

-0.16987266E  00 

j 

i 

2 

7 

O.OOOOOOOOE  00 

0.28B02817E  OX 

0.1888 1666E  00 

i 

2 

XO 

O.OOOOOOOOE  00 

-0.28S02398E  01 

O.OOOOOOOOE  00 

'j 

3 

8 

-0.89601S81E  OX 

O.OOOOOOOOE  00 

0.16987266E  00 

3 

8 

O.OOOOOOOOE  00 

0.60093021E  00 

0.16X05571E  00 

j 

3 

XX 

O.OOOOOOOOE  00 

-0.60092365E  00 

-0.500000X2E  00 

LOADING 

CONDITION  3 

HUB 

VX 

VY 

6 

X 

0.10882807E-08 

-0.9X6X6930E-05 

0 .5 3 78862 8 E— 03 

1 

2 

0.6105176 OE— 05 

— 0. 18308965E— 08 

0.1078X268E-02 

3 

O.OOOOOOOOE  00 

-0.2X38XXX8E-08 

0.125 08582 E-02 

HUB  ELEMENT 

MX 

MY 

S 

X 

X 0.X500X078E  03 

O.OOOOOOOOE 

00 

-0.36885769E 

02 

1 

5 -0.15001829E  03 

O.OOOOOOOOE 

00 

0.39612683E 

02 

1 

6 O.OOOOOOOOE  00 

0.61621985E 

03 

0.26673182E 

02 

X 

9 O.OOOOOOOOE  00 

-0.6 162  1 8B  7E 

03 

— 0.250000C3E 

02 

2 

2 — 0.68315930E  03 

O.OOOOOOOOE 

00 

-0 .2 1 1 1626  IE 

02 

2 

3 0.68315112c  03 

O.OOOOOOOOE 

00 

-0.27736138C 

02 

2 

7 O.OOOOOOOOE  00 

0. 1639675  IE 

06 

0.73868686E 

02 

2 

XO  O.OOOOOOOOE  00 

0.X6396731E 

06 

— 0.25000003E 

02 

3 

8 — 0.58880173E  03 

O.OOOOOOOOE 

00 

-0.32265872E 

02 

3 

8 O.OOOOOOOOE  00 

0.17061191E 

06 

0.89992 126E 

02 

3 

XX  O.OOOOOOOOE  00 

- 

0.1706 116 9E 

06 

-0.25500003E 

02 

ELEMENT  SIGMA  0 

SIGMA  8 

SIGMA  MAX 

SIGMA  P 

X 

0. X01639E  OX 

0.636262E 

OX 

0.605137E 

01 

O.OOOOOOE 

00 

X 

-O.X63X76E-01 

-0.805666E-02 

O.OOOOOOE 

00 

O.OOOOOOE 

00 

X 

0. 100007E  01 

0.655636E 

OX 

0.606060E 

ox 

O.OOOOOOE 

00 

3 

0.896239E  OX 

0.339627E 

ox 

0.709878E 

ox 

O.OOOOOOE 

00 

3 

-0.a0989XE-02 

0.397365E-01 

O.OOOOOOE 

00 

O.OOOOOOE 

00 

3 

0.895833E  OX 

0. 365599E 

01 

0.712197E 

ox 

O.OOOOOOE 

00 

5 

0. X01681E  OX 

-0.267939E 

02 

O.OOOOOOE 

00 

O.OOOOOOE 

00 

9 
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00 
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00 

9 
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02 

O.OOOOOOE 

00 
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00 

6 
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02 
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00 

6 
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a 
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02 
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9 
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02 
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02 
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02 
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02 
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00 
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00 
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00 
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02 
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00 
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92 

Results  for  Stress  Analysis 
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RUCKLING  load  COEFFICIENT  = 0.6275E  Ct> 


buckling  mode  shape 


THETAX 
C.1633E-01 
-0.1625E-01 
“0*251 1 E-03 
0.1653E-01 
-0.1645E-01 


THETAY 
-0.1645E-01 
-0.1625E-01 
-0.2512E-03 
0.1653E-01 
0 • 1633E-0 1 


0 • 2 5 0 9 F 
• 2 A 3 1 E 
O.ICOOE 
0.2536E 
G.25G9F 


T mF TAX 
-C.1173E-G3 
^ • 3 2 A I E-G 1 
-G.3225E-CI 
-0.1173E-C3 


thftay 
-0.3225E-C1 
-0.  il73E-02 
“0 • 1 »7AE“C3 
0.32A1E-C1 


Lowest  buckling  mode  of  grillage 


UNLIMITED 


General  Buckling  of  the  grillage  with  only  the  first  5 modes  shown 


BUCKLING  LOAD  COEFFICIENT 


BUCKLING  MO OF  SHAPE 

THETAX  TriETA 

-0.5310E-04  0.5310E 

0.5310E-U4  G.5310F 

C.1440E-39  0.3765E 

-0.5310F-04  -C.5310F 
0.531CE-Q4  — 0 » 53  1 OF 


0.8 12br-C3 
■0 .8  126E— 03 
U.3239E-C2 
0.312bF-C3 
0.8 I26E-C3 


BUCKLING  l.OAL)  COEFFICIENT 


BUCKLING  MODE  SHAPE 

THETAX  THETAY 

-0.3935F-04  0.5025E-04 

0.8901E-Q5  -C.3249E-04 
0.1042E-03  -0.6077E-G4 
0.890 1 E— 05  -0.3249E-04 
-0.3935E-04  C.5025E-04 


TnETAY 
0 . 1 9 k 9E  — 
0.2786L-C'* 
0.27e6E-04 
0. 1909E-0H 


U.8819F-03 
0.2321E— 03 
0.5468F-G9 
0 .232  IE-03 
U.8820E-03 


BUCKLING  LOAD  COEFFICIENT 


BUCKLING  MODE  SHAPE 

THETAX  THETAY 

0 . 32  33E-04  0.9099E— 05 

— 0.5021E— 04  -0 • 3948E-04 
0.6118E-0A  0 • 1039E-03 

-0.5021E-0A  -0.39A8E-0A 
0.3233E-0A  0.9099E-05 


THETAX 
0.2805E— OA 
■0.1922  E-OA 
■C.1922E-CA 
0.2805E— OA 


0. 1239L 
0. 7296L 
0.7296L 
O'.  1239L 


BUCKLING  LOAD  COEFFICIENT 


BUCKLING  vOD 
THETAX 
0.2505E-0A 
0.250AE-0A 
— 0.9565E— 10 
-0.2505E-0A 
-0.250AE-0A 


SHAPE 

THETAY 

■O.250AE-0A 

O.250AE-0A 

0.A2UE-09 

0.250AE-0A 

0.250AE-0A 


ThLTAY 
0.3561E-10 
0. 7760E—0A 
0. 7759L-CA 
0.21&3E-09 


THETAX 

0.7760E-GA 

G.6973E-1C 

0.13H7E-09 

0.7759E-04 


0. 39 lOt 
0.339  ft 
0.K2A8E 
0.6521c. 


BUCKLING  LOAD  COEFFICIENT 


BUCKLING  MODE  SHAPE 

THETAX  THETAY 

-0.6553E-05  0.7871E-05 

-0.2517E-OA  0.2555E-0A 

-0.2235E-04  0.2355E-0A 

-0.2517E-04  0.2555E-CA 

-0.655AE-05  0.7870E-05 


W THETAX 

0 . 909  3E-03  0.A309E-OA 

0.2  36  IE— CA  0.7862E-05 
0.9060E-08  0.7863E-05 

0.2361E-0A  0.A309E-0A 

0.9093E-03 


0.2607E— 03 
0 . 2h7AE— 03 
0-2A75E-03 
j . 2 bOb  L— 0 J 
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APPENDIX  E 


RUNNING  TIMES 
■*.  IBM  1130 

a.  Stress  Analysis:-  3 hubs,  9 beams  - 5 mins  20  secs 

b.  Lowest  Buckling:-  9 hubs,  24  beams  - 37  mins  46  secs 

c.  General  Buckling:-  9 hubs,  24  beams  - 15  mins 

d.  Natural  Frequencies:-  9 hubs,  24  beams  - 6 mins  34  secs 

2.  UNI  VAC  1108 


a. 

b. 

c. 

d. 


Stress  Analysis:- 
Lowest  Buckling :- 
General  Buckling:- 
Natural  Frequencies :- 


3 hubs,  9 beams  - 1 CP  sec 
9 hubs,  24  beams  - 1.5  CP  secs 

9 hubs,  24  beams  - 2.5  CP  secs 

10  hubs,  11  beams  - 1.3  CP  secs 


3.  KPF-9 

a.  Stress  Analysis  (in-core ):- 

b.  Lowest.  Buckling :- 

c.  Stress  Analysis  (out  of  core):- 

d.  Data  Generation:- 

e . Data  Plot : - 


4 hubs  - 38  secs 
9 hubs  - 23  secs 
30  hubs  - 2 min  24  secs 
12  hubs  - 20  secs 

18  hubs  - running  time  - 2 min  23  sec 
plotting  time-  11  min 
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APPENDIX  F 

RESTRICTIONS  ON  THE  SIZE  OF  THE  PROBLEM  AND  RESTRICTIONS  ON  THE 


ELEMENTS  THAT  CAN  BE  USED 


Maximum  number  of  elements  = 80 
Maximum  number  of  loading  conditions 


Stress  Analysis  (in-core  solution) 

NBEAM  (21+17NL0AD)  + 1 8NL0AD  + 3 JOINT  ( 3 JOINT+NLOAD ) < 1 600 

Lowest  buckling  mode 

20  NBEAM  + 2(3J0INT)2  <1600 

Stress  Analysis  (out  of  core  solution) 

NBEAM ( 56+1 9NL0AD)  + 6NL0AD  + 9J0INT  + 2IBAND  + 

3MAXIB  (9MAXIB  + NLOAD)  <1600 
and 

3 JOINT  (3 JOINT  + NLOAD )<  200000  , 

Data  Generation 

26NBEAM  + 7NBEAM  NLOAD  + NLOAD  + 2J0INT  <16000 
Data  Plot 

7NBEAM  (3  + NLOAD)  <16000 


UNLIMITED 


Restrictions  on  the  Elements  that  can  be  Used 


Element 
■ Type 


Facility 


Stress  Analysis 
Lowest  Buckling  Mode 
General  Buckling  Mode 
Natural  Frequencies 


APPENDIX  & 


Vo, wo 

Ur  j 


MITED 


Nil m 


rn«tt 


nsjH  ♦ | 


Lowest  buckliwot 

ROUTINE  CkCHl) 

rfStfH  * O 

6UCKLTH  Gr  Koo  TZt4  e 


°*<i  r*M  mj 

man  **i  U<* 


c.H  fKoMU 

A«  i^iat 


call  e«  rr 


wss  fcooTlrft 


C*LL  C«  IT 


■1 

■ 

Cm*| 

4i  cnUoMt 
<4mu*s 

NOTATION 


APPENDIX  H 


End  restraint  constants  for  single  ended  in  plane  element 

<*,  - -1  v tb  * -*>’**■  sb  * -°,wb> 


*b 

*d 

•f 

*hb,  hf,  hw 
m 

P 

Vx,  Vy 
W 

**,  y 

*xi * yi 
**2 9 y 2 


width  of  faying  flange 
depth  of  frame  web 
width  of  flange 

thickness  of  faying  flange,  flange  and  ftrame  web  respectively 

A concentrated  load  acts  at  a distance  mB  from  end  0 of  the 
element 

Uniform  load  applied  to  the  element,  positive  downwards 
Rotations  in  x and  y directions 

Displacement  perpendicular  to  the  plane  of  the  grillage, 
positive  downwards 

Rectangular  coordinate  axes  at  each  hub  point 
Distances  from  hub  to  end  0 of  element 

Distances  from  hub  to  end  B of  element 


*Z 


A 

A1,  B1 

*B 

E 

G 

I 

J 

Mx,  My 

NBEAM 

NT 

*N1 

N2 

JOINT 

S 

T 

W 

0 

SIGMA  0, 
SIGMA  M 
SIGMA  P 


Distance  of  the  outermost  fibre  of  the  beam  cross-section 
from  the  neutral  axis 

Effective  shear  area  of  element 

Cross  sectional  area  and  length  of  pillar  support 
Length  of  an  element 
Young's  Modulus 
Modulus  of  shear 

Moment  of  inertia  of  cross-section  for  bending  out  of  the  plain 
Torsional  constant 

Moment  in  x and  y directions  respectively 
Total  number  of  elements  in  grillage 
Parameter  to  determine  version  of  the  program 
Label  number  of  element  (at  end  0) 

Element  type  number 
Number  of  hubs 

Foroe  perpendicular  to  the  plane  of  the  grillage 
Torsional  moment 

Concentrated  load  applied  to  the  element,  positive  downwards 
Angle  element  makes  with  x axis 
Angle  between  element  and  image 
SIGMA  B Stress  at  end  0 (or  end  B)  of  element 

Stress  at  a maximum  in  bending  moment  distribution  if  any 
Stress  at  point  of  application  of  concentrated  load  if  any. 


Notation  marked  with  • is  illustrated  in  figure  6 
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